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Abstract 

In this work a transformation acoustics scheme for generic elastic media is developed. Our 
approach starts form the decomposition of the elasticity tensor in terms of its eigentensors, 
an idea previously used by Norris. While Norris' transformation acoustics is restricted to the 
special class of so-called pentamode materials, we show that a similar scheme can be defined 
for the most general elasticity tensor. As in case of Norris' model (and in sharp contrast 
to transformation optics), the compatibility equations of the transformation medium are not 
purely algebraic and it is not guaranteed that solutions to these equations exist for any choice 
of material parameters and coordinate transformation. Nonetheless, it is shown that our 
scheme yields new cloaking solutions for certain classes of materials. In particular, we present 
the first application of a transformation based device for a non-scalar wave equation outside 
of the field of electromagnetics. 

1 Introduction 

Transformation optics [IH3] is one of the most important design tools in the field of artificial elec- 
tromagnetic materials (metamaterials) . A unique combination of characteristics makes this tool 
extremely powerful: the trajectories of light in the metamaterial are simply defined as geometric 
deformations (coordinate transformations) of the trajectories in free space; the solutions of the 
Maxwell equations in the metamaterial together with its material properties are obtained from the 
free space solutions and the deformations in a completely algebraic way; once certain intuitive re- 
strictions on the deformations are respected, all interfaces between free space and the metamatrial 
are automatically refiectionless [I]. 

This combination for the first time allowed to derive the necessary material parameters of 
almost arbitrary optical devices in a simple, intuitive way. Not surprisingly, besides the original 
proposal of a cloak [TJ[1] a pletora of new concepts was proposed (see, for example, [31[5HI]). 

Of course, these successes raised the question, whether similar tools can be defined for other 
theories of physics, notably different types of waves. The first attempt to apply similar concepts 
to acoustic waves [5] appeared in the same year as the electromagnetic cloak and became an in- 
dependent field of research [9lfT2] . Still, the manipulation of acoustic waves turns out to be much 
harder than of its electromagnetic counterpart. The reason for this lies in a very specific charac- 
teristic of the Maxwell equations, which is not shared by most other wave problems in physics: the 
existence of a premetric formulation. Indeed, by a suitable choice of electromagnetic variables, the 
complete Maxwell equations can be written in a form that does not include any spacetime metric 
or any object derived from the spacetime metric [13] • Therefore, in these variables the Maxwell 
equations — and in consequence any solution thereof — do not carry any information about the 
local structure of spacetime or the local choice of coordinate^. This peculiar feature of electrody- 
namics makes transformation optics possible: starting with a free space solution of the Maxwell 
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1 This fact is often referred to as form invariance of the Maxwell equations. This term can be misleading, since 
— in a suitably covariant notation — any sensible law of physics is form invariant under diffeomorphisms. 
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equations in coordinates x and with spacetime metric g, one applies a coordinate transformation 
to obtain the same solution in coordinates x with spacetime metric g; since the Maxwell equations 
do not depend on the metric, one can easily interpret the transformed solution as a solution in the 
original coordinates x with spacetime metric g, which is exactly the transformation optics solution 
that belongs to the coordinate transformation x —> x(x). Though mostly used in the context of 
deformations of the free space solutions, transformation optics can be applied to any constitutive 
law, including nonlinear |14U15) and nonlocal (16) media relations. 

The theory of acoustic waves in elastic bodies is not premetric and thus the scheme sketched 
above cannot apply to this theory in general. Still, it was found [51[T7] that in certain special cases 
transformation acoustics exists. Consider a limit of the full theory that can be simplified to a 
scalar field <j>{x,t) obeying a Helmholtz equation 

M = 4>- (i) 

In a covariant notation this equation reads 

DrfWrf = (2) 

where D is the covariant derivative with respect to the spatial metric 7. The second covariant 
derivative just acts on a scalar and thus simply can be replaced by a plain derivative d. The first 
covariant derivative constitutes the gradient of a vector 

Drf , v l = ^d j( j> , (3) 

which can be written as 

£y = ^=MV70 > (4) 

with 7 = det(7ij). Therefore, the Helmholtz equation becomes 

-j=di (y/rftdrf) = m'uKg = '<j>. (5) 

This equation appears to be metric independent if k and M lJ can be interpreted in a suitable way 
as material parameters. It is important to realize that k does not transform as a scalar, but as a 
relative scalar and, similarly, M is a tensor density rather than a tensor. To our knowledge, all 
transformational design tools outside of transformation optics known so far are based on Eq. ([5]). 
Consequently, all these design tools are restricted to scalar wave problems. It is the purpose of 
this paper to suggest a route to overcome this restriction, whereby linear elastodynamics will serve 
as a guiding example. 

The idea presented here is a generalization of Norris' pentamode transformation acoustics [10 . 
In this model, Norris' defines a special class of metamaterials in elastodynamics, whose wave 
equation essentailly reduces to Eqs. ([T])-©. After a brief review of linear elastodynamics (Sec. [2]) 
and the special case of pentamode materials (Sec. [3]) we will show in Sec. 2] that the wave equation 
of a large class of elastic media can be cast into the form 

^ l\ 1 IhC-'j'i /•>,/'./ Pi, (6) 

J 

where the indices /, J take values (1,2, ... , 6). This means, that the tensorial wave equation of 
elastodynamics in these cases can be rewritten as a system of six coupled scalar wave equations. In 
consequence, the transformation acoustics scheme based on Eq. ([S]) conceptually can be applied to 
this large class of media straightforwardly, but the ensuing compatibility equations for the material 
parameters turn out to be complicated and — as already found by Norris [10) — non-algebraic. 
The concept of Eq. ([6]) can be further generalized to include all linear elastic media by adding 
appropriate source terms to that equation. Despite their complexity it is shown in Sec. [5] that new 
solutions for linear radial transformations in cylindrical and spherical coordinates can be found 
based on our concept, which includes new types of cloak solutions. 
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2 Linear elastodynamics 



Our model are the linearized elasticity equations [THHH] in terms of the stress tensor cr = (cr 4 - 7 ), 
the infinitesimal strain tensor e = (ey) and the elasticity tensor C = (C l ° kl ). In order to fix our 
notation we provide a brief overview of the relevant equations in this section. 

The stress tensor is subject to momentum and angular momentum balance. Due to the latter, 
cr must be symmetric, <7 V} — a 3 " 1 , while the former can be cast into the form 

D j( r ij + P F i = P W i . (7) 

Here, p{x) is the density of the volume element, pF l are the total volume forces and W % is the 
acceleration of the volume element. In this paper we always assume F = 0. Starting from Eq. 
([7]) leads to a very restricted version of transformation acoustics, in particular, the special case 
of inertial transformation acoustics [TTK^nidl] is not covered at all. To get more flexibility, a 
anisotropic mass density is indispensable (see Ref. [TU] and references therein). Therefore, for the 
purpose of transformation acoustics Eq. (JT]) is replaced by 

PM J • (8) 

The infinitesimal (or small) strain tensor is the Lie derivative along the displacement vector u 
of the spatial metric, 



2^.7 , dj = 2 i D i u 3 + D J u i) > ( 9 ) 



where D is the covariant derivative with respect to the spatial metric 7. The elasticity tensor 
defines a linear constitutive law between the stress and the strain tensor: 

o*j = C^ kl e M (10) 

Finally, the acceleration vector W must be related to the displacement vector field by 



d 2 u 

W = ^=u. (11) 



Due to the characteristics of the stress and strain tensor C must be symmetric in its first two 
and in its last two indices. Thus, C has 36 independent parameters and may be represented as 
a generic 6x6 matrix. Additional constraints on the elasticity parameters emerge if an elastic 
potential W with 

a y = — , C ljkl = , (12 

detj deijdeki 

shall be defined. In that case C %]kl = C kl%3 and in analogy to electromagnetism [T3] such a material 
can be called skewonless. As already mentioned by Milton and Cherkaev [35], the elasticity tensor 
of any skewonless material can be decomposed at|j 

6 

c m = ^ KiS ii s u (13) 

7=1 

This decomposition is not unique. The scalars Kj can be chosen as eigenvalues to the normalized 
eigentensors S l - ' , Tr(S ■ S) = 1, alternatively they can be fixed as numbers K] = ±1, 0. In general, 
the Ki are scalar fields. In contrast to Ref. [32] the decomposition is not fixed in this article. 



2 In contrast to indices of physical space, no summation convention is invoked for the space of eigentensors and 
eigenvectors, labeled by capital Latin indices. 
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3 Norris' transformation acoustics revisited 



An important precursor to results derived in this work is the paper by Norris jlOj . Thus, it appears 
helpful to reformulate the results of that paper in our notation in a first step. A more careful 
re-derivation following the logic and notation of this paper can be found in Ref. [23 . 

Norris considers a very special material, a so-called pentamode medium. In terms of the decom- 
positon (I13[) this class of materials is characterized by the fact that it has only one nonvanishing 
eigenvalue: 

K\ = K K 2 = K 3 = K i = K 5 =K 6 = (14) 

This is also referred to as a material with five "easy modes" , since five out of six modes of the 
stress tensor are not related to the modes of the elasticity tensor. 

In addition, Norris requires that the eigentensor S which corresponds to the nonvanishing 
eigenvalue K must have vanishing divergence: 

DiS ij = (15) 

As shown by Norris, the wave equation simplifies in a very interesting way for this special system. 
A quantity with dimension of a pressure can be defined by contracting S with the strain tensor e: 

p = /\.S"<, ,, = -KS lj D lUj . (16) 

With this definition, the stress tensor becomes o % i = —S^p and, consequently, the differential 
equation of the stress tensor (O can be written in the compact form 

= -S ij D jP = // ,/•;•' . (17) 

The time derivative of the displacement field can be related to the time derivative of p by differ- 
entiating Eq. (|16l) with respect to time 

p = -KS ij DiUj , (18) 

where we assumed that the elasticity tensor is not time dependent. Combining (fTTf and ([T8"]) one 
arrives at a wave equation for p: 



P 



KDi (s l ip^S kl ) Dip = l\ !),(•"> I), /> (19) 



It is seen that for the special case of a pentamode medium the wave equation of elasticity theory 
reduces to a scalar wave equation of the "pseudo-pressure" p. The tensor G = S ■ p 1 ■ S can be 
interpreted as effective metric, but it is important to realize that the covariant derivatives are still 
defined with respect to the spatial metric -y rather than G. Turning the covariant derivatives into 
plain derivatives (cf. Eqs. ([I])-©) yields 

P = ^di^jG^djp . (20) 
V7 

Consider a coordinate transformation x =>■ x(x) applied on the wave equation (j 191) - Since we 
use a covariant formulation the effect of the transformation boils down to putting a bar over each 
non-scalar quantity: 

p = KDi {S ik p^S lj ) D 3 p = KDiG ij Djp , (21) 
which in the language of Eq. (|2"0"|) becomes 

p = ^=di^G ij d jP . (22) 
V7 
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To conclude the program of transformation acoustics, the transformed wave equation should be 
re-interprctcd in the original coordinates, which at the level of Eq. ((22|) means that the barred 
partial derivatives should be replaced partial derivatives without a bar. However, to obtain the 
material parameters of the transformation acoustics device one has to get back to a wave equation 
in a covariant formulation, 

Equating Eq. (|2"Tj) or (f^^j) with Eq. (|23p unveils the necessary material parameters of the trans- 
formation acoustics device: 



K=^K S lJ p^S kl = ^S l 3p-^S kl (24) 



At this point, one has to remember the constraint (|15[) . which played a crucial role in the derivation 
of the wave equation (|19l) . Indeed, after the coordinate transformation S satisfies the constraint 
DiS lJ = 0, but in a transformation acoustic interpretation the altered material matrix has to 
satisfy DiS l i = 0. Therefore, the second equation in ((24)) requires an appropriate redefinition 
of the elasticity tensor and the mass density matrix of the transformation acoustics material. As 
discussed in Ref. [TU] this in general can only be achieved if an anisotropic mass density is assumed. 

Though the above equations in principle define transformation acoustics for any pentamodc 
material and any coordinate transformation, an important restriction has to be stressed: In con- 
trast to transformation optics, the compatibility equations for the material parameters of the 
transformation medium are not entirely algebraic, but include the differential equation DiS 1 -* = 0. 
The formalism of transformation acoustics does not guarantee that a solution to this equation can 
be found. 



4 Transformation acoustics for generic media 

It is the aim of this section to show that a transformation scheme similar to Norris' theory of 
transformation acoustics conceptually can be defined for a much larger class of materials. In 
particular, we relax the constraint (|14l) and in principle allow for an arbitrary set of eigenvalues 
Kj. 

4.1 The general divergence free material 

In a first step we derive a transformation acoustics scheme that follows closely the idea by Norris. 
Consequently, we assume that for a certain decomposition (1131) all Si obey 

DiSy =0 V/ . (25) 

Obviously, this constrains the possible spatial dependence of the elasticity tensor in a quite drastic 
way. In Cartesian coordinates one simply concludes DiC l ^ kl = diC^ kl = ^{diKi)S) 3 S kl . 

The derivation of the wave equation closely follows the one of the previous section. For each 
matrix Si a "pressure" can be defined as 

pi = -KiS l /D lUj . (26) 

This allows to rewrite the stress tensor as a lJ = — YIiPi^i 3 and thus the momentum conservation 
becomes 

D j o*t=-52S?D jPl = p i j & . (27) 
i 
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This equation can be solved for ii, the result is used in Eq. (|2l))) twice differentiated with respect 
to time: 



Pi 



= K I sr i D kPm ]D l Y,sy PJ 

j 

K l D k Sf n p£&iD j p J (28) 



= KiDi&ijDtfj 
j 

As a result, the wave equation of elastodynamics has been transformed into six coupled wave 
equations of six scalars pi. If the mass density matrix p is assumed as symmetric, these coupled 
wave equations are characterized by 21 effective metrics Gu = Gji. 

Despite the rather complicated form of the wave equations ([28)) , this reformulation allows to 
proceed the transformation acoustic scheme along the lines of Norris' pentamode transformation 
acoustics. Applying the coordinate transformation x — > x(x), Eq. (f2"5f becomes 

pj = £ I\ 1 1)/,"/ j I) ,i>.i = J2 ^diV7G%d jPj . (29) 

To reinterpret the result as mimicking of the coordinate transformation, the barred partial deriva- 
tives are again changed back to partial derivatives without bar. To obtain the correct covariant 
derivatives in the coordinate system x, the known trick (|24l) of rescaling by the square root of the 
metrics is used: 

P 1 = £ ^d^GyjdjPj 



Ki^D k y-LsrP^SyD jPj (30) 



£" 

J 



kxDi (s* m p-\syD 3 pj 



Thanks to the reformulation of the wave equation in terms of the scalars pi, the rescaling method 
used by Norris (Eq. (|23])1 also defines a transformation acoustics scheme for a generic divergence 
free material. However, the material parameters of the transformation medium are not defined in 
a unique way, but the possible realizations of the transformation medium are defined by the set 
of equations 

6 

Qijkl = ^K I S i I j Sf , (31) 
i=i 

Kj = ^Kj , (32) 

~Sfp- k l~S l i = ^~Sfp~ k rS l ] , (33) 

DiSf =0 VI . (34) 

Obviously, there exist solutions to the two algebraic constraints (|3"2"|) and (f3"51) for any choice of 
material parameters and coordinate transformation. However, as in the case of Norris' pentamode 
transformation acoustics, the last constraint is not an algebraic equation and it is not guaranteed 
that a solution exists. Still, a generic transformation acoustic scheme cannot be expected, since 
elastodynamics, in contrasts to electrodynamics, is not a premetric theory. Therefore, it is seen 
as a success of the reformulation in terms of the scalars pi that the consequences of the missing 
premetricity of the theory can be cast into the simple form of Eq. (IMt . As will be shown below, 
the recipe derived in this section can be helpful in certain situations to obtain transformation 
acoustic results for materials different than just pentamode media. 
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4.2 Decomposition of the general material 

Let's simply drop the condition that the matrices Si should be divergence free. It was seen in 
the previous section that a reformulation of the wave equation in terms of the scalars pi was the 
important step to bring the tcnsorial wave equation of clastodynamics into a form that is accessible 
to the standard techniques of transformation acoustics. Since the definition of pi in Eq. (|26D is 
purely algebraic, these scalars exist completely independently of the constraint (|25[) . However, the 
simple form of the momentum conservation as given in Eq. (|27p is no longer correct; instead, that 
equation has to be replaced by 

=-5>, (S?pi) = - J2(j} Pl + S?D 3PI ) = //,,/' , (35) 



where we introduced the source terms DiSj 1 — ji. From the definition of the scalars pi in Eq. 
(l2l)l) one obtains 



pi = —KiSfDiUj = KiS l /D t (^J k x Dt {Sfpj^j 
= £ Ki ( Di S? - jj) P~l (SfDt + j k j) PJ . 



(36) 



As is seen, the reformulation in terms of the scalars pi still allows to reformulate the wave equation 
of clastodynamics in terms of six coupled differential equations in the six scalars pi. While in the 
case of the divergence free material, these coupled equations where all of Helmholtz type, the 
equations (|36p include additional inhomogeneous terms, which make the situation more complex. 
Despite this fact, one can proceed straightforwardly in the program of transformation acoustics. 
Since all quantities of Eq. (|3"6"|) are just scalars, vectors or tensors the effect of a coordinate 
transformation x — 5- x is merely to put a bar over each quantity. Since all covariant derivatives 
hit either a scalar or constitute the divergence of a vector, this barred equation again can be 
rewritten as 

* = E Ki (7^^ - & + PJ ■ (37) 

The transformation acoustic scheme following from this transformation again interprets the barred 
partial derivatives as derivatives without bar. Transforming back to covariant derivatives without 
bar, the transformation acoustic analogue to Eq. (|37[) is 



Pi 



r r (38) 

J2^Ki{D t Sy-ji)^(SfD l+ ^ Pj . 



Comparison of the Eq. (|38|) with Eq. (1361) unravels that the material parameters of an acoustic 
transformation material mimicking the coordinate transformation x — > x follow from the set of 
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equations 
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Qijkl = ^K I S i I i Sf , (39) 
1=1 

Ki = ^=K l , (40) 
V7 

S}%lS l j = ^_Sfp- k yS l ] , (41) 

jitf® = ^j k r P - k rs l J , (42) 

^Pilfj = ^LfrPilfj , (43) 

US'/ fi ■ (44) 



The new inhomogeneous terms yield two new algebraic conditions in Eqs. (|42l) and (|43[) . Still, 
solutions to all algebraic conditions (f4Tfl) - (H5)) can be found for any choice of material parameters 
and coordinate transformation. More importantly, the inhomogeneity enters the differential def- 
inition of jj in Eq. (|44[) . which makes it very difficult in general to find a valid solution to all 
equations (|to" ]t -([i3 |) . 

An interesting simplification of the above equations occurs, if all Si are either full rank or 
divergence free. This constitutes a restriction on the choice of material parameters, which however 
is less restrictive than the divergence free materials discussed in the previous section. In this 
particular case one can define the alternative source terms 

ju:= !° fc . if AS? = o, 

1 SYijD^S/ otherwise, 

and the wave equations (|3^|) and (f3"T| become 

pj = J2Ki(D l - JIt )S?p H 1 S l ](D 3 +jj j )pj 
j 

= Y, K i^ d ^-Ji*)S} k p k - l 1 S l j(d J +jj 3 )pj (46) 
J V ' 

where the effective metrics Gu are defined as in Eq. (f^5|) . The last equation differs from the 
respective form of the wave equation of the divergence free case only by the replacement d — > 
d ± ji. Evidently, the source terms do not affect the transformation acoustic interpretation of 
Ki and Gij. Therefore, the material conditions of the acoustic metamaterial, which corresponds 
to the coordinate transformation x — > x, follow in analogy to the divergence free case, Eqs. 
(35M33), as 

6 

Qijkl = ^K I S) j Sf , (47) 
i=i 

Ki = , (48) 

V7 

~Sfp- k yS l j = ^Sfp^S 1 ] , (49) 

£>>"■' H" - ifA5??=0 ' (50) 
otherwise. 
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Since this formulation of material equations suppresses any dependence of the algebraic conditions 
on the source terms ji , it provides considerably more freedom to find solutions to these conditions 
that obey the definition of jj = ji i from Eq. ([501) . 



5 Radial transformations and the cloak 

In many applications the transformation device is either of spherical or cylindrical shape and one 
is interested in modifying the trajectories of waves exclusively in the radial direction. The best 
known examples are the spherical [T] and cylindrical [24] cloak. Thus, we want to investigate in 
this section if and how purely radial transformations can be extended by means of the formalism 
developed in this work. Various strategies for acoustic cloaks have been presented in the literature 
(see Ref. [12] and references therein), but to our knowledge all these proposals assume an idealized 
situation where a single scalar wave is governed by a Helmholtz type equation. Since our formalism 
should be able to cover systems with more than one scalar mode it should be possible to find 
entirely new cloaking applications. 



5.1 Radial transformations in cylindrical coordinates 

Let us start by analyzing the situation of a cylindrical system. In cylindrical coordinates the 
metric tensor is given by 

= diag(l,r 2 , 1) , tf*] = diag(l, 1 1) . (51) 

All information about the local structure of space and the choice of coordinates is encoded in this 
object. An important quantity in our formalism is the covariant divergence of the matrices Si, 
which is obtained from 

DiSj = diSj + T l mi Sj ] + Y 3 mi SY n . (52) 
In cylindrical coordinates, the Christoffel symbols 

Fli = \l" n (dilmk + d klml - d mlkl ) (53) 

are given by 



diag(0,-r,0) , T e = U , T z = , (54) 




which yields the divergence (remember that Si must be symmetric) 

DiSf = diSf + -SV - rS e r e , (55) 
r 

DiSf = diSf + -Sf , (56) 
r 

D t S l i z = diSf + -SF . (57) 
r 

As has been worked out in the previous sections, the tour de force of our transformation 
acoustics scheme is to relate the divergence of Si after the coordinate transformation, DiSf , to 
the divergence of a suitable choice of new material parameters Si, DiSf . Therefore, we develop 
in a first step a strategy to solve this problem for radial transformations in cylindrical coordinates. 
Consider a linear radial transformation applied to our system, 

f = Ar + B, 6 = 6 , z = z . (58) 
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Under such a transformation, the metric and the Christoffcl symbols transform as 

[%} = diag(-^, (r ~f )2 , 1) , = diag(A 2 , ^ 1) , (59) 

^ 0\ 

diag(0,-(f-i?),0) , ^=1^ r , f* = 0. (60) 

0/ 

Therefore, the divergence of the matrices Si is given by 

DiSf = BiSf + ^-5^ - (f - S)£f , (61) 

r — B 

D t S? = 9i5f + — ^— , (62) 
r — B 

DiSf = BiSf + r-^SY . (63) 
r — B 

To solve the constraints (l4"4")l or ([5"0"f new matrices Si must be defined in such a way that 

DiS^ can be related to DiS 13 . Ideally, DiS lJ is chosen componentwise proportional to DiS l] : 

D t S ij = aWiS ij no sum over j. (64) 

In general, this will lead to a complicated system of differential equations. In the example at 
hand, however, a simple solution can be found if we assume that all Si are diagonal in cylindrical 
coordinates. This means that we have at most three (linearly independent) matrices Si. The three 
independent Si correspond to the three pressure modes, while the three shear modes become "easy 
modes" . Since this material is shear free and in general has an anisotropic mass density as well as 
an anisotropic elasticity tensor, it can be seen as a generalization of the anisotropic metafluid [25 , 
for which a transformation acoustics scheme is known already. 

With this simplifying assumption the covariant derivatives (|61[) - (|63l) become 

DiSY = d r SV + —^—SY - (f - B)S 6 i 6 , (65) 
r — B 

DiSf = BeSf , (66) 

DiSf = 3Jr . (67) 

The simple form of the covariant derivatives allows to define a linear transformation, Si — T ■ S/, 
such that the divergence with respect to D is transformed into a divergence with respect to D. 
Indeed, the structure of the dependence in S rr suggests the ansatz 

Si = exp[/(r)]Sj , (68) 

where /(r) follows from 

d r f(r) = ^—-- =* /(r)=ln— -. (69) 
r — B r r 

An additional scaling factor is needed to correct the S® 9 term in (1551 and thus we find from the 
rescalings 

- = r_-B_ B „ ^ - ee = (r~B) 2 §ee § = r-B_ B (7Q) 

the new covariant derivatives 

Dr grr = r _£ Dr srr ^ Dg §ee = (r - B) ^ D j„ = r-B_ D ^ § „ (n) 
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The rescaling of the Sf z components is not fixed by Eqs. (|65 p — (|67 |l but has been chosen for 
convenience. 

Since all matrices Si are diagonal, the above rescalings indeed may written as the matrix 
equations 

r-B (r- B) 2 r - B 
S I = T-S I = S I -T , T = diag( , { - ^-L, ). (72) 

Most importantly, the matrix T has been constructed in such a way that 

D i T i j S jk = T k jDiS ij . (73) 

The construction performed so far solves the differential constraint (|44l) . From this result, the 
remaining constraints Eqs. (|40j) - (|43[) are easily solved. Indeed, from the "commutation relation" 
(l73l) it follows that jj in Eq. (|44|) is obtained by the same rescaling matrix as S^: 

jj = T'i'j) (74) 

This, in turn, implies that on the left hand sides of Eqs. PT1)- P3"|) the same resaclings are applied 
and thus these three equations are simultaneously solved by 

p=^T-p-T. (75) 
V7 

Eqs. (j?2"]) and ([75)1 define the material parameters of a transformation acoustics scheme for a 
coordinate transformation of type (|58p and for a material with three eigentensors Si. These 
tensors have to be diagonal in cylindrical coordinates but need not be of full rank or divergence 
free. 

5.2 Radial transformations in spherical coordinates 

In spherical coordinates (r, 9, </>) the metric tensor reads 

[ 7<J -] = diag(l, r 2 sin 2 <j>, r 2 ) , [ 7 «] = diag(l, —^-r, \) , (76) 

r z sm </> r z 

which yields the Christoffel symbols 

V din-iO. -rsin 2 (p,-r) , T e = [ ± 6 cot^] , T z = ( - sin cos ^ 0| .(77) 

Thus, the divergence of the matrices Si in these coordinates reads 

DiSf = diSf + -SY + cot (j)Sf r - r sin 2 4>Sf - rSf^ , (78) 
D t Sf = diSf + -S} 6 + 3 cot 0Sf , (79) 
= fltSf + ^ + cot - sin cos «£Sf 9 . (80) 
In analogy to the previous section we apply a linear radial transformation 

f = Ar + B, 9 = 9, <$> = <$>• (81) 

The metric tensor under this transformation changes to 

m = diag(i,, ^ sin 2 0, ^) , [^] = diag(A 2 , {f _ £ ^ ^ J^) , (82) 
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As in cylindrical coordinates, all factors A 2 drop out in the calculation of the Christoffcl symbols, 
such that they are obtained from the expressions in Eq. (|77p simply by replacing r everywhere by 
f — B. Consequently, the divergence of Si in spherical coordinates is given by 

DiS} r = diSf + —t—SY + cot 0Sf r -(f-B) sin 2 $Sf - (f - B)Sf , (83) 
r — B 

Di Sf = B t Sf + Sf + 3 cot 4>Sf , (84) 

DSf = BiSf + y^Sf + cot 4>S^ - sin 4> cos $S$ e . (85) 

In order to reinterpret this result in the original coordinate system x the same strategy as in the 
previous case proves successful: First, we demand that all Si are diagonal in spherical coordinates, 
which simplifies Eas. (jg3"|) - (l8"5]) to 

DiSf = B r SY + yz^SY -(f-B) sin 2 - (f - B)^ , (86) 

DiSf = B e S d i 9 , (87) 
DiSf = B^Sf + cot $Sf - sin $ cos 0Sf e . (88) 

Again, a diagonal matrix can be defined, which relates between the divergence with respect to 7 
and with respect to 7. Instead of the rescaling matrix T in Eq. ((72|) . a slightly different matrix U 

r~B (r — B) 2 (r - B) 2 
S 1 = U-S I = S I -U , [/ = d iag(— — , 1 r2 > , [ r2 > ) . (89) 

must be used here, which also fulfills the important "communtation relation" (|T3")) . Therefore the 
material parameters in case of a spherical transformation are obtained from Eq. (|89H and 

r, 1 • p = ^u-p-u. (90) 

V7 



6 Summary and conclusions 

In this work we developed a generic theory of transformation acoustics applicable to any elasticity 
tensor. In contrast to electromagnetics, elastodynamics is not a premetric theory and thus it is not 
expected that a transformation based design scheme exists for any choice of material parameters 
and coordinate transformations. Our approach to attack this problem is based on the previous 
works by Milton et al. [S] and Norris [TU] . Following these authors we reformulate the wave equation 
in terms of the moduli Ki and the eigentensors Si of the elasticity tensor. In this formulation 
the complex tensorial wave equation turns into six coupled wave equations in six scalar quantities 
Pi. In consequence, for each of these equations a transformation scheme similar to the one of the 
Helmholtz equation can be invoked. As a result, compatibility equations for the eigentensors and 
moduli are found according to which the transformation material has to be defined. 

Of course, the metric dependence of the theory is preserved by this reformulation and thus the 
related problems with transformation acoustics cannot be eliminated. In our formulation, these 
limitations result in non-algebraic compatibility equations (see Eqs. (| 3 2 1) - (|3"4"1) . Eqs. (|4TH (|4"4"|) and 
Eqs. (|48 p — (|50 p ). Still, as a main result of this article we find that thanks to our reformulation 
the complicated dependence of the tensorial wave equation on the spatial metric can be hidden 
completely and that the metric dependence can be reduced to simple conditions on the divergence 
of the six eigentensors Si. To our knowledge this scheme for the first time provides a realistic 
transformation tool outside of electromagnetics that is applicable to more than just a single scalar 
mode. 

As an application of our algorithm we studied linear radial transformations in cylindrical and 
spherical coordinates, which include the cylindrical and spherical cloak. It was shown that in 
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this special case transformation media can be found for any material with three pressure modes 
but vanishing shear modes (anisotropic liquid). This result constitutes a generalization of previous 
studies of liquids in transformation acoustics, where the anisotropy of the metafluid only concerned 
its mass density tensor |25) . 

Our result suggests that the transformational design approach could be a powerful tool in a 
much broader context than known today. The basis for such applications seems to be a clever re- 
formulation of the wave equation in such a way that the required coordinate transformation easily 
can be implemented. A similar result has recently been developed by Garcia-Meca et al. [26]: these 
authors show, though for a different system than in this article, that with a suitable choice of vari- 
ables certain transformation acoustics schemes become possible or exhibit different solutions than 
known before. Still, the results found so far may be insufficient in various applications. While in 
electromagnetics transformation devices can be matched without any reflections to the surround- 
ing space A\ , interface conditions between standard elastic media and acoustic metamaterials from 
transformation acoustics have not been studied extensively so far. 
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